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SIMULATION CAPABILITY FOR DYNAMICS OF TWO-BODY 
FLEXIBLE SATELLITES* 


Fred Austin** and George Zetkov 
Grumman Aerospace Corporation 
Bethpage, New York 


Abstract 


An analysis and computer program were prepared to 
realistically simulate the dynamic behavior of a class of 
satellites consisting of two end bodies separated by a 
connecting structure. The shape and mass distribution of 
the flexible end bodies are arbitrary; the connecting struc- 
ture is flexible but massless and is capable of deployment 
and retraction. Fluid flowing in a piping system and rigid 
moving masses* representing a cargo elevator or crew 
members, have been modeled. Connecting Structure 
characteristics, control systems, and externally applied 
loads are modeled in easily replaced subroutines. Sub- 
routines currently available include a telescopic beam- 
type connecting structure as well as attitude, deployment, 
spin, and wobble control. In addition, a unique mass bal- 
ance control system was developed to sense and balance 
mass shifts due to the motion of a cargo elevator. The 
mass of the cargo may vary through a large range. Numer- 
ical results are discussed for various types of runs. 

I. Introduction 

Often, when a new satellite configuration requires in- 
vestigation, the equations of motion must be derived, pro- 
grammed, and checked in order to predict the dynamic 
behavior of the vehicle. In this paper, a new computer 
program is described which eliminates these tasks for 
a class of two-body flexible satellites. One application 
is the rotating counterweight space station. 

The end bodies are referred to herein as the Labor a - 
tory and the Counterweight , whereas the entire struc- 
ture is referred to as the Space Station. The structure 
separating these bodies is referred to as the Connecting 
Structure . Both the Laboratory and the Counterweight 
are modeled as flexible structures with arbitrary shape 
and mass distribution (see Figure 1). The Connecting 
Structure is also flexible, but its mass is neglected. 



Fig. 1 Idealization and Coordinate Systems 


Laboratory, and the effect of fluid pumped through 
a piping system on the Laboratory. An example 
of a realistic configuration which may be treated 
is shown in Figure 2. Various features including 
control systems and Connecting Structure char- 
acteristics have been included in subroutine form 
so that these items may be easily replaced with 
minimum disruption to the main program. A 
subroutine is also provided for programming 
forces and torques; these may be applied at any 
mass point on the Space Station. Special con- 
straint options permit the user to rigidize the 
entire Space Station or certain portions of the 
vehicle; thus the easier-understood rigid-body 
results can be compared with the more complex 
flexible-body solutions. 


Problems which may be studied using the program in- 
clude deployment and retraction of the Connecting Struc- 
ture with simultaneous or sequential spin-up and spin- 
down, the effects of moving rigid masses such as a 
cargo elevator or crew members on board the 


It is also possible to use the program to study 
the dynamics of the Laboratory alone when no 
Counterweight is present. Thus, many rotating 
and nonrotating satellite configurations can be 
studied since the Laboratory characteristics may 
be varied as input data. 


*This work was performed for the NASA Langley Research Center under contract NAS 1-10973. The complete study 
is documented in References 1 and 2. Most of the content of this paper was presented at the AIAA Dynamics 
Specialist Conference at Williamsburg, Virginia on March 1973 (AIAA Paper No. 73-320, entitled "Simulation 
Capability for Dynamics of Rotating Counterweight Space Stations"). 


** Structural Mechanics Engineer 
^Guidance and Control Engineer 


The authors appreciate the technical and administrative assistance provided by Dr. R. Fralich, the NASA 
Project Monitor. Mr. S. Goldenberg performed all theoretical work required for the preprogram, which 
synthesizes the modes of the main structures. Dr. J. Markowitz's many significant technical contributions to 
this program are appreciated. Mr. E. Lowe programmed and checked both the main program and the preprogram, 
and his excellent work is gratefully acknowledged. Many of Mr. J. Smedfjeld's useful comments were incorporated 
into the manuscript. 
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Fig. 2 Typical Configuration That Can Ec Analyzed Using the Computer Program 


IL Idealization and Method of Analysis* 

The coordinates used in the analysis are shown in 
Figure 1. The X axes are a system of mean axes moving 
with the average or ’’rigid-body’’ motion of the Laboratory 
and are used as the reference system for its elastic mo- 
tion. Similarly, the Y axes are a set of mean axes 
for the Counterweight. ( R } and the Euler angles 
{ 7 } locate the X axes in space, and ( R } and j v \ 
locate the Y axes with respect to the X axes. The 
elastic linear and angular deflections with respect 
to the mean axes are j and {# jj, respectively, 
for a typical mass point m\ (i = 1, . . . , n) on the 
Laboratory and{q a )and ( 0 a [for n typical mass point 

m a (a - 1, n) on the Counterweight. These 

elastic coordinates were linearized whenever there 
was an analytical or computational advantage to do 
so; however, the rigid -body coordinates | R | 7 j, 

{ H j t and { v } are not linearized so that large rigid- 
body motions may be studied. Instead of Euler - 
angle rates ({7} and {*) ^angular velocities (jtu } 
and{<t|YJ) of the X and Y axes, respectively, are used 
as integration coordinates in order to simplify the 
equations of motion. The angular-velocity com- 
ponents are not derivatives of generalized coordinates 
but are derivatives of quasi-coordinates. (3) 

Newton and Euler equations of motion were 
written for an arbitrary number of lumped masses 
on the Laboratory and the Counterweight**, and the 
position and orientation of these masses is also 
arbitrary. The rotatory inertia of each structural 
mass point is included, as is linear and angular 
momentum of a fluid confined within a pipe segment 
on a typical Laboratory mass point, Any Labora- 
tory mass points may also include the two fluid- 
system reservoirs; one is nominally filling and 
the other is nominally emptying. Uniform 
flow is assumed; thus slosh is not allowed. 

A maximum of eight moving rigid masses (a cargo 


elevator, crew motion, etc. ) arc idealized as 
point masses with no rotary inertia. Their 
rigid body motion is a prespecified function of 
time; however, the additional motion due to 
structural deformation of the path is estimated 
by averaging the motion of the surrounding 
structural masses and in this way is approx- 
imately included in the analysis. 

The equations of motion were first written in vec- 
tor form, and then converted to a convenient matrix 
form by using a set of identities developed for this 
purpose (see Appendix B of Reference 1). Next, 
modal or constraint relationships (or a combina- 
tion of the two) were used to reduce the number of ■ 
coordinates. For example, one type of constraint 
that may be selected by the program user is to ri- 
gidize the Laboratory or the Counterweight. After 
reduction of coordinates, the equations are linearly 
combined into a reduced number of equations by an 
automated technique described in Appendix D 
of Reference 1. The result is identical to that 
w'hich would be obtained by writing Langrange *s 
Equations for a system containing quasi -coordin- 
ates; however, the technique used in this analysis 
is more direct and simpler. The equations now 
have a symmetric acceleration coefficient matrix 
which was used to facilitate computation as 
follows. Because of nonlinear effects this ma- 
trix, and consequently its inverse, varies; there- 
fore, the equations must be solved at each time 
point. Since more rapid computing algorithms are 
available for solving equations with symmetric 
coefficient matrices, the symmetric form led to 
appreciable time savings. Also, less storage is 
required since only the upper triangle of the matrix 
must be stored. 


♦Details of the analysis are presented in Reference 1. 

**In the computer program the maximum number of lumped masses is 100 for each body. 
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When modes are used to represent the elastic 
motion of the Laboratory or Counterweight, a 
number of free-free elastic modes (but not rigid- 
body modes) of that structure are read into the 
program. Omitting the rigid-body modes uniquely 
determines the location of the X and Y mean axes 
of Figure 1 (see Section 4. 3, 5 of Reference 1). 

The modal masses and frequencies are required 
by the program and are used to define the stiff- 
ness characteristics of the structure; thus it is not 
required to read in a large physical stiffness matrix. 
Modal damping is assumed, and a different percent- 
age of critical damping for each mode may be used. 

HI* Obtaining the Modes of Vibration 

A separate preprogram was prepared to synthe- 
size the modes of the Laboratory or Counterweight 
from a knowledge of the modes of their various com- 
ponent substructures. The synthesized modes may 
be passed directly from the preprogram to the main 
program. 

To provide the versatility for synthesizing a 
wide variety of configurations, the procedure is de- 
veloped for the general seventeen- module idealiza- 
tion shown in Figure 3. The user has the capability 
of eliminating any of the modules from tills most 
general arrangement so long as he does not discon- 
nect the structure. In this way, many configurations 
having a lesser number of flexible modules may be 
studied. The general configuration in Figure 3 allows 
up to five modules to 1x3 appended to one another in an 
arbitrary manner, and up to ten modules to be appen- 
ded to one core module. 
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Fig. 3 General 17 -Module Configuration for Modal 
Synthesis 


One feature of the synthesis procedure is that the sub- 
structure modal matrices maybe supplied in coordin- 
ate systems that are not parallel to the coordinate system 
in which the results, the coupled modes, are obtained. 
Accordingly, modules need not be in the same plane. In 
fact, they may be skewed at any angle in space. 

Another feature of the procedure permits the user to 
supply constrained substructure modes. These are modes 
which were obtained for idealizations where constraints 
were employed; for example, in a beam analysis axial 
extension may have been neglected. 


Capability is not provided to automatically synthesize 
structures where appendages axe interconnected; however, 
modes for these structures may be obtained by other pro- 
cedures and then used as input data to the main program. 
Thus, the main program is not restricted by the limita- 
tions of the preprogram. 

IV. Control Systems 

Five major types of control systems have been de- 
veloped and programmed in subroutine form. The 
controt and associated command subroutines were 
originally developed on a rigid-vehicle idealization 
and were then implemented on the flexible -vehicle 
computer program. Controt actuators and sensors 
were assumed to move with the total motion of the mass 
point to which they were attached. The control 
systems arc described briefly below. 

Three-Axis Attitude -Position Control fwhen the vehicle 
is not spinning) . Sixteen jets are mounted on the Labora- 
tory for attitude-position control. The attitude error e 
and error rate c are measured about each axis, and 
the e and e combination is used to decide which jets 
to fire. 

Spin Rate Control . The jets are also used for spin- up 
and spin- down maneuvers when these maneuvers are 
commanded by a spin command subroutine. In addi- 
tion, these jets can be automatically fired to maintain 
a constant spin rate within a certain tolerance 
established by the control-system dead band. 

Counterweight Position Control. This control system 
is used for deployment and retraction maneuvers. Un- 
like the other control systems, the actuator dynamics 
are not simulated. Rattier, an idealized motor is 
assumed which is able to precisely control the unde- 
formed length 4$(t) of the connecting structure; i.e., 
is set equal to a pre specified function of time 
contained in a position command subroutine. The in- 
put data to this subroutine are: the initial and final 
positions, the time for deployment to begin, the mag- 
nitude of an on-off constant acceleration f X | , and a 
maximum velocity,/ . u 

0 max 

Wobble Control. A control moment gyroscope (CMG) 
is used to damp undesirable wobble due to gyroscopic 
effects when the vehicle is rotating. The highly ef- 
ficient 90° h-lag law(4~7) is used to accomplish this 
task. 

Mass Balancing for Spinning Vehicle . When the cargo 
elevator is in motion, the shifting of a large mass causes 
undesirable motions of the spinning vehicle. To balance 
tills effect, an accelerometer detects the resulting 
accelerations, and a balance mass is moved to cor- 
rect both the center of mass shift and the cross pro- 
ducts of inertia of the vehicle. 

Miscellaneous Command Subroutines . The same sub- 
routine used to command the undeformed length of the 
Connecting Structure is also used to command the posi- 
tion (along each of the three axes) of every moving 
point mass except the balance mass (which is governed 
by the balance mass control system). There are a 
maximum of seven such moving masses which may be 
used to simulate the cargo elevator, crew members, 
etc. Another subroutine is used to prespecify the 
fluid velocity in the piping system, i.e. , an idealized 
pump is assumed. Automatic shutoff occurs when a 
reservoir is either emptied or filled. This subroutine 
also generates the spin-rate command used as a 
reference input by the spin- rate control system. 
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V. Miscellaneous Program Capabilties 

For reference, the program computes the position 
of the Space Station center of mass in inertial coor- 
dinates, the total angular momentum vector projected 
onto inertial coordinates, and the total system kinetic 
energy. The output includes J 6 )and ( * } , the linear and 
angular displacements, respectively, of the Connecting 
Structure attachment point on the Laboratory relative 
to the attachment point on the Counterweight. In addi- 
tion, the program can compute the internal resultant 
force and torque vectors on any surface that separates 
the structure Into two free bodies. This is accomplished 
by the acceleration method, which appears to be more 
accurate than the stiffness-damping matrix method 
for the case where a truncated number of modes repre- 
sents the deformation of the vehicle. 

VI. Program Checkout 

Several test problems were run to check the pro- 
gram. In the test configuration, the Laboratory was 
idealized using eight structural masses and the Coun- 
terweight idealization had five structural masses. A 
telescopic boom was used as the Connecting Structure, 
First the entire vehicle was commanded to be rigid, 
and the time history was found to be identical to results 
obtained using the w r ell-known analytical solutions for 
a rotating rigid body. Next, a separate program was 
written by an independent programmer using the theo- 
retical expressions of Reference 1. However, this 
check program was much simpler than the main pro- 
gram since the expressions were applied only to the 
test configuration. Also, the time- and storage -saving 
manipulations used in the main program were not used 
in the test program. Several test runs were made, in- 
cluding vibration excited by: giving each variable and 
its derivative a different initial condition while rotating; 
fluid motion; point masses in motion; and a retraction 
maneuver. Time histories of the two programs agreed 
in each case. 

The preprogram which synthesizes the modes of the 
Space Station was also tested by comparing results for 
several configurations with independently obtained 
results. The independent results were obtained by com- 
puting a mass and stiffness matrix for the entire struc- 
ture by a standard technique and applying an existing 
eigenvalue program. Agreement was achieved in each 
case. 

VIL Numerical Results Demonstrating 
the Computer Programs 

Configuration . The NASA Langley Research Center 
provided Grumman with the configuration shown in 
Figure 2 for the purpose of demonstrating the com- 
puter program. The Laboratory is composed of a 
central core module with two relatively rigid appen- 
dages and two solar panels which arc very flexible. 

The Connecting Structure is a telescopic beam which 
can be fully retracted. The Counterweight is com- 
posed of three relatively rigid modules. The mass 
of the entire vehicle is approximately 77, 600 kg, and 
its deployed overall length is approximately 78.4 m. 

For additional descriptive details of this configuration 
see Reference 1. 

First, the modes of the core and appended modules 
were computed by Grumman, based upon stiffness 
and mass properties generated by North American 
Rockwell. The modes of the solar panels were 
computed by Fairchild Industries and Wolf Research. ' 
These modes were used as input data to the pre- 


program in order to obtain the free-free modes of the 
Laboratory and the Counterweight. The Counterweight 
was relatively rigid. The lowest flexible Counterweight 
frequency was 6.851 Hz, whereas the sixth Laboratory 
frequency was .1382 Hz, It was therefore decided to 
idealize the Counterweight as a rigid body in the time- 
history computer program. The Laboratory was 
idealized using 72 mass points. Six flexible laboratory 
modes were then used to reduce the number of 
numerical- integration coordinates. Because of the 
relatively high flexibility of the solar panels, most of 
the motion in these modes is solar-panel motion. 

A total of 18 coordinates are used in the time-history 
runs. These are the six laboratory modes, six 
rigid-body coordinates locating the mean axes for the 
Laboratory, and six rigid-body coordinates for the 
Counterweight. 

The connecting structure was assumed to be an 
elastic tubular beam with uniform characteristics 
per unit length. To approximate telescopic charac- 
teristics of the beam during deployment and re- 
traction, new stiffness properties are computed at 
each time interval based on the beam length. 

All control system jets, the CMG, and the control- 
system sensor package were mounted on the Labora- 
tory core module. 

A detailed description of the configuration is 
presented in Section G. 0 of Reference 1. 

Selection of Runs. The runs which will be discussed 
arc attitude control, deployment, spin-up, wobble 
control, elevator motion with mass balancing, and 
fluid being pumped between reservoirs. It is also 
possible to use the program to study many of these 
effects occurring simultaneously; however, this was 
not done in the present study. Since the runs 
performed were selected primarily to demonstrate 
the capability of the computer program, there was 
no attempt to optimize any of the operation param- 
eters such as control-system gains. Also, in the 
case of the deployment and spin-up maneuvers, the 
jet thrusts were increased to unrealistic values in 
order to complete the maneuvers within 45 sec, 
thereby saving computer time. 

Although the runs described herein are primarily 
flexible-vehicle runs, in each case the rigid-vehicle 
run was also performed. Rigid izati on was accom- 
plished by using the program's constraint option. 

Attitude Control. During the attitude-control 
maneuver, the system was not rotating and the 
Connecting Structure was fully retracted. The 
attitude control system developed for the rigid 
vehicle required no modification for use on the 
flexible-vehicle idealization. The three components 
of the Euler angles {7} orienting the vehicle are 
shown in Figure 4, Curves for a rigid-body run 
were overlayed with the flexible body curves of 
Figure 4 and no difference could be discerned. 

The Space Station was initially tilted so that each 
Euler angle was .01745 rad (1.0 deg). The control 
system then reduces each angle to the commanded 
value of zero. Similar behavior occurs along each of 
the three axes. The jets first apply a torque to begin 
correcting the attitude angle. Then a torque is 
applied in the negative direction to slow down the 
Space Station's angular rate. Linear deformation 
in the Connecting Structure is shown in Figure 5. 

Since the Connecting Structure is fully retracted, 
the illustrated deformation actually respesents the 
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During Attitude Control Maneuver 


very small deflections of the relatively stiff Con- 
necting Structure docking hatch. This connection 
is only approximately represented as discussed in 
Section 4.4. 2.1 of Reference 1. 


Deployment with Spin-Rate H old . Figure 6 illus- 
trates a deployment maneuver at .2 rpm. All motion 
in this run is in the spin plane. Figure 6A shows the 
prespecified motion of the undcformed Connecting 
Structure. The Connecting Structure deploys from 
an undeformed length l of 0 to 42. 822 m. If the 
spin control system were not operational, the spin 
rate would decrease to maintain constant momentum; 
however, the command to maintain a constant spin 
rate was given during this run. To accomplish the 
deployment within 40 sec, the jet thrust was increased 
from a nominal value of 222. 4 N (50 lb) to 4448. 2 N 
(1000 lb). Figure 6B shows the axial deformation in 
the beam during deployment. Throughout the maneu- 
ver, the slow spin rate causes a slight centrifugal 
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Fig. 6 Deployment Maneuver at . 2 RPM 


force. Therefore, there is a small tendency for the 
end bodies to separate naturally; however, this would 
occur at a very slow rate. The deployment command 
is given at t - 1 sec, and the Laboratory and the 
Counterweight are pushed apart (£03 >0 as shown in 
Figure 6 A) causing compression in the beam. When 
the maximum velocity £ 03 = 1.270 m/sec (50 in. /sec) 
is reached, and deployment proceeds at constant 
velocity, the beam is expanded slightly by the centrif- 
ugal force. At approximately 35 sec, the deceleration 
begins ( t 03 < 0 ) and the expansion in the beam is 
increased significantly. Deployment ends at about 
40 sec. The final expansion is much larger than 
the initial compression, mainly because the Con- 
necting Structure is more flexible when more of it 
is deployed. For the same reason, the transient 
vibration occurs at a lower frequency when deploy- 
ment is near completion. Figure 6 C illustrates 
the bending of the Connecting Structure during deploy- 
ment. This bending occurs partially as a result of 
the Coriolis forces, but primarily it is due to the spin 
jets torquing the Laboratory to maintain a constant 
angular velocity. 


Spin-Up. Figures 7 and 8 illustrate a spin-up 
maneuver. - The Space Station is initially rotating at 
. 2 rpm and, at 5 sec, the command is given to in- 
crease the spin speed to 4 rpm. All motion in this 
run is in the spin plane. In order to accomplish the 
maneuver in 40 sec, the jets on the Laboratory were 
increased from 222.4 N (50 lb) to 66.7 kN (15, 000 lb). 
Figure 7 A illustrates the increase of the spin speed, 
and Figure 7B illustrates the corresponding increase 
in the axial extension of the Connecting Structure due 
primarily to the centrifugal force. Figure 7C shows 
the bending in the Connecting Structure during the 
spin-up maneuver. Unusually high bending deforma- 
tions occur as a result of the Large torques on the 
Laboratory generated by the increased jet thrusts. 
Figure 8 shows the largest component of the deforma- 
tion at mass point 17, a very flexible point on a solar 
panel, and the Largest component of the force exerted 
by the solar panel on the core module at the root of 
the panel (see Figure 2). 

Quiescent State. When the Space Station is rotating 
in its nominal state of pure spin (i. e. , in a state of 
pure rotation about the Xj axis with no vibration), 
constant elastic deformations occur due to the cen- 
trifugal force. This state is known as the quiescent 
state. During the runs which were made when the 
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Fig. 7 Spin Speed and Cantilever Deflection of 

Connecting Structure During Sptn-Up Maneuver 
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Fig. 8 Selected Deflection and Load Components 
During Spin-Up Maneuver 


Space Station was rotating at its nominal spin speed 
of 4 rpm, the initial conditions are a variation from 
the quiescent state* Before making these runs* the 
quiescent state deformations were determined by 
setting all of the damping coefficients (for both the 
Laboratory and the Connecting Structure) to 80% of 
their critical values. A short run was made, and the 
deformations rapidly damped to their quiescent 
values. The quiescent deformations were highest 
at certain points on the solar panels. As an example, 
the deformation in the X 3 direction at mass point 17 
(q 17j 3 ) was approximately 540 mm (21.4 in.). 


Wobble Control . Figure 9 illustrates the perfor- 
mance of the wobble control system. Initially, the 
deformations were set to their quiescent values, and 
the second component of was given a wobble 
component of .001 rad/sec. Up to approximately 
27 sec the curves are essentially identical to a run 
made for a rigid Space Station. This indicates the 
usefulness of the mean axes; one reason that they 
were used was that they move at the average motion 
of the deformed system. After 27 sec, some small 
higher-frequency oscillations predominate due to 
elastic vibration. The only modification required to 
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Fig. 9 Components of Angular Velocity Illustrating 
Wobble Control System 


the CMG control system was to increase the amount 
of wobble at which the CMG stops trying to control 
wobble. Before this modification was made, the CMG 
sensor reacted to the residual vibration and the CMG 
continued to operate in the wobble-damping mode 
throughout the entire run. 


Elevator Motion with Balance Mass Control . In 
the runs described in this paragraph, a 4,530 kg 
(10,000 lb) elevator and a 2,270 kg (5, 000 lb) balance 
mass are initially located on the X 3 axis, at X 3 ---19.4 
and 8.9 m, respectively. The elevator moves towards 
the balance mass at t=5 seconds. All rigid-body and 
flexible motion is in the spin plane. Initially, the 
Space Station is rotating in the quiescent state. The 
first balance-mass control system which was devel- 
oped operated properly on a rigid-vehicle idealiza- 
tion; however, using the current computer program 
it was found that vehicle flexibility coupled undesir- 
ably with the control system. During this run the 
vehicle vibrated at a high-frequency vibration 
associated with the Connecting Structure axial mode, 
and the overall response of the balance mass was 
sluggish and unstable. The satisfactory control 
illustrated in Figure 10 was achieved after 
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considerable modification of the control law. The 
control system is required to balance an elevator 
motion of 7.00 m. The curve showing the error in the 
position of the Space Station center of mass indicates 
that there is a lag in the response of the balance mass; 
however, by the end of the run the Space Station is 
balanced. The deformations in the Connecting 
Structure during this run are shown in Figure 11. 

Note that there is an initial vibration although the 
elevator does not begin its motion until 5 sec. This 
initial vibration results because the initial quiescent 
deformations were obtained from a run where the 
elevator and balance masses were not present. 

The main Connecting-Structure bending effects are 
caused by the Coriolis forces exerted by the moving 
masses on the Laboratory and by the spin jets 
which torque the Laboratory to maintain the com- 
manded constant spin speed. The deformation at 
mass point 17, which is on a solar panel, &nd the 
internal force exerted by the solar panel at its root 
arc shown in Figures 12A and B. 


in Figure 13. Rimping proceeds until this reservoir is 
empty at t = 35. 7 see, when the pump suddenly shuts 
down. After shutdown, fluid remains in the pipeline. 

All motion during this run occurs in the spin plane. 
Figure 14 shows the deformations of the Connecting 
Structure. The bending which occurs during pumping 
is illustrated in Figure 15. The primary reason for 
this bending is that the resultant of the Coriolis forces 
exerted by the fluid acts to the left of the cm of the 
Laboratory as illustrated in the figure. This force also 
tends to slow the spin speed of the Laboratory slightly 
(from .4189 to .4149 rad/sec.). The control system 
is not present during this run; therefore, the spin 
jets do not turn on to correct the spin speed. 
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Fig. 13 Fluid Motion on the Laboratory 



Fig, 11 Cantilever Deflection of Connecting Structure 
During Mass Balancing 
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Fig. 14 Cantilever Deflection of Connecting Structure 
During Fluid Motion 
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Fluid Pumped Through the Laboratory. In this 
run only, the Laboratory contained a fluid system con- 
sisting of a pipeline connecting two reservoirs. Unlike 
the system illustrated in Figure 2, the centerlines of 
both reservoirs and the pipe were located on the 
axis. Pumping begins at 2 sec. The fluid velocity and 
height of the fluid in the .emptying reservoir are shown 
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Fig, 15 Sketch Illustrating Bending (Exaggerated) in 
Connecting Structure During Fluid Motion 
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Abstract 

An analysis and computer program were prepared to 
realistically simulate the dynamic behavior of a class of 
satellites consisting of two end bodies separated by a 
connecting structure. The shape and mass distribution of 
the flexible end bodies are arbitrary; the connecting struc- 
ture is flexible but massless and is capable of deployment 
and retraction. Fluid flowing in a piping system and rigid 
moving masses, representing a cargo elevator or crew 
members, have been modeled. Connecting Structure 
characteristics, control systems, and externally applied 
loads. are modeled in easily replaced subroutines. Sub- 
routines currently available include a telescopic beam- 
type connecting structure as well as attitude, deployment, 
spin, and wobble control. In addition, a unique mass bal- 
ance control system was developed to sense and balance 
mass shifts due to the motion of a cargo elevator. The 
mass of the cargo may vary through a large range. Numer- 
ical results are discussed for various types of runs. 

1. ’ Introduction 

Often, when a new satellite configuration requires in- 
vestigation, the equations of motion must be derived, pro- 
grammed, and checked in order to predict the dynamic 
behavior of the vehicle. In this paper, a new computer 
program is described which eliminates these tasks For 
a class of two-body flexible satellites. One application 
is the rotating counterweight space station. 

The end bodies are referred to herein as the Labora- 
tory and the Counterweight , whereas the entire struc- 
ture is referred to as the Space Station. 1 he structure 
separating these bodies is referred to as the Connecting 
Structure . Both the Laboratory and the Counterweight 
are modeled as flexible structures with arbitrary shape 
and mass distribution (see Figure 1). 1 he Connecting 

Structure is also flexible, but its mass is neglected. 

Problems which may be studied using the program in- 
clude deployment and retraction of the Connecting Struc- 
ture with simultaneous or sequential spin-up and spin- 
down, the effects of moving rigid masses such as a 



Fig. 1 Idealization and Coordinate Systems 


Laboratory, and the effect of fluid pumped through 
a piping system on the Laboratory. An example 
of a realistic configuration which may be treated 
Is shown in Figure 2. Various features including 
control systems and Connecting Structure char- 
acteristics have been included in subroutine form 
so that these items may be easily replaced with 
minimum disruption to the main program. A 
subroutine is also provided for programming 
forces and torques; these may be applied at any 
mass point on the Space Station. Special con- 
straint options permit the user to rigidize the 
entire Space Station or certain portions of the 
vehicle; thus the easier-understood rigid-body 
results can be compared with the more complex 
flexible-body solutions. 

It is also possible to use the program to study 
the dynamics of the Laboratory alone when no 
Counterweight is present. Thus, many rotating 
and nonrotating satellite configurations can be 
studied since the Labor atory characteristics may 
be varied as input data. 


cargo elevator or crew members on board the 
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Fig. 2 Typical Configuration That Can Be Analyzed Using the Computer Program 


II. Idealization and Method of Analysis* 

The coordinates used* in the analysis are shown in 
Figure 1. The X axes are a system of mean axes moving 
with the average or "rigid-body" motion of the Laboratory 
and are used as the reference system for its clastic mo- 
tion. Similarly, the Y axes are a set of mean axes 
for the Counterweight. | R \ and the ICider angles 
j 7 ) locate the X axes in space, and | R } and { V f 
locate the Y axes with respect to the X axes. The 
elastic linear and angular deflections with respect 
to the mean axes are | cjj j and j 9 j), respectively, 

for a typical mass point (i ” 1 n) on the 

Laboratory and |q a )nnd | Tor a typical mass point 

m a (a - 1, n) on the Counterweight. These 

elastic coordinates were linearized whenever there 
was an analytical or computational advantage to do 
so; however, the rigid-body coordinates | R } , { 7 ), 

{ R | , and j * 1 are not linearized so that large rigid- 
body motions ma}' be studied. Instead of Eulcr^ 
angle rates ( (y \ and j ) ^angular velocities ({w ' | 
and ( W Y)) 0 f the X and Y axes, respectively, are used 
as integration coordinates in order to simplify the 
equations of motion. The angular-velocity com- 
ponents are not derivatives of generalized coordinates 
but are derivatives of quasi-coordinates. (*0 

Newton and Euler equations of motion were 
written for an arbitrary number of lumped masses 
on the Laboratory and the Counterweight**, and the 
position and orientation of these masses is also 
arbitrary. The rotatory inertia of each structural 
mass point Is included, as is linear and angular 
momentum of a fluid confined within a pipe segment 
on a typical Laboratory mass point. Any Labora- 
tory mass points may also include the two fluid - 
system reservoirs; one is nominally filling and 
the other is nominally emptying. Uniform 
flow is assumed; thus slosh is not allowed. 


elevator, crew motion, etc. ) are idealized as 
point masses with no rotary inertia. Their 
rigid body motion is a prespecified Function of 
lime; however, the additional motion due to 
structural deformation of the path Is estimated 
by averaging the motion of the surrounding 
structural masses and in this way is approx- 
imately included in the analysis. 

The equations of motion were first written in vec« 
tor form, and then converted to a convenient matrix 
form by using a set of identities developed for this 
purpose (sec Appendix B of Reference 1). Next, 
modal or constraint relationships (or a combina- 
tion of the two) were used to reduce the number of 
coordinates. For example, one type of constraint 
that may be selected by the program user is to ri- 
gidizc the Laboratory or the Counterweight. After 
reduction of coordinates, the equations are linearly 
combined into a reduced number of equations by an 
automated technique described in Appendix D 
of Reference 1, The result is identical to that 
which would be obtained by writing Langrange's 
Equations for a system containing quasi-coordin- 
ates; however, the technique used in this analysis 
is more direct and simpler. Hie equations now 
have a symmetric acceleration coefficient matrix 
which was used to facilitate computation as 
follows. Because of nonlinear effects this ma- 
trix, and consequently its inverse, varies; there- 
fore, the equations must be solved at each time 
point. Since more rapid computing algorithms are 
available for solving equations with symmetric 
coefficient matrices, the symmetric form led to 
appreciable time savings. Also, less storage is 
required since only the upper triangle of the matrix 
must be stored. 


A maximum of eight moving rigid masses (a cargo 
♦Details of the analysis are presented in Reference 1. 

♦♦In the computer program the maximum number of lumped masses is 100 for each body. 



When modes are used to represent the elastic 
motion of the Laboratory or Counte rweight, a 
number of free-free elastic modes (but not rigid - 
body modes) of that structure are read into the 
program. Omitting the rigid-body modes uniquely 
determines the location of the X and Y mean axes 
of Figure 1 (see Section 4.3,5 of Reference t). 

Hie modal masses and frequencies are required 
by the program and are used to define the stiff- 
ness characteristics of the structure; thus it is not 
required to read in a large physical stiffness matrix. 
Modal damping is assumed, and a different percent- 
age of critical damping for each mode may be used. 

EL Obtaining the Modes of Vibration 

A separate preprogram was prepared to synthe- 
size the modes of the Laboratory or Counte rweight 
from a knowledge of the modes of their various com- 
ponent substructures. The synthesized modes may 
be passed directly from the preprogram to the main 
program. 

To provide the versatility for synthesizing a 
wide variety of configurations, the procedure is de- 
veloped for the general seventeen-module idealiza- 
tion shown in Figure 3, Hie user has the capability 
of eliminating any of the modules from this most 
general arrangement so long as he does not discon- 
nect the structure. In this way, many configurations 
having a lesser number of flexible modules may be 
studied. The general configuration in Figure 3 allows 
up to five modules to be appended to one another in an 
arbitrary manner, and up to ten modules to be appen- 
ded to one core module. * 
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Fig. 3 General 17 -Module Configuration for Modal 
Synthesis 


One feature of the synthesis procedure is that the sub- 
structure modal matrices may be supplied in coordin- 
ate systems that are not parallel to the coordinate system 
in which the results, the coupled modes, are obtained. 
Accordingly, modules need not be in the same plane. In 
fact, they may be skewed at any angle in space . 

Another feature of the procedure permits toe user to 
supply constrained substructure modes. These are modes 
which were obtained for idealizations where constraints 
were employed; for example, in a beam analysis axial 
extension may have been neglected. 


Capability is not provided to automatically synthesize 
structures where appendages are interconnected; however 
modes for these structures may be obtained by other pro- 
cedures and then used as input data to the main program. 
Thus, the main program 16 hot restricted by the limita- 
tions of the preprogram. 

IV. Control Systems 

Five major types of control systems have been de- 
veloped and programmed In subroutine form. The 
control and associated command subroutines were 
originally developed on a rigid-vehicle idealization 
and were then implemented on the flexible -vehicle 
computer program. Control actuators and sensors 
were assumed to move with the total motion of the mass 
point to which they were attached. The control 
systems are described briefly below. 

Hi re e -Axis Attitude -Position Control (when the vehicle 
is not spinning) . Sixteen jets are mounted on the Labora- 
tory for attitude -position control. The attitude error e 
and error rate e are measured about each axis, and 
the e and c combination is used to decide which jets 
to fire. 

Spin Rato Control . The jets are also used for spin-up 
and spin-down maneuvers when these maneuvers are 
commanded by a spin command subroutine. In addi- 
tion, these jets can be automatically fired to maintain 
a constant spin rate within a certain tolerance 
. established by the control- system dead band. 

Counterweight Position Control. Hiis control system 
is used. for deployment and retraction maneuvers. Un- 
like the other control systems, the actuator dynamics 
are not simulated. Rather, an idealized motor is 
assumed which is able to precisely control the unde- 
formed length ^ 0 (t) of the connecting structure; Le. , 

^o(t) Is set equal to a prespecified function of time 
contained in a position command subroutine. The in- 
put data to this subroutine are: the initial and final 
positions, the time for deployment to begin, the mag- 
nitude of an on-off constant acceleration | }' [ , and a 
maximum velocity^ . u 

0 max 

Wobble Control. A control moment gyroscope (CMG) 
is used to damp undesirable wobble due to gyroscopic 
effects when the vehicle is rotating. The highly ef- 
ficient 00° h-lag law(4-7) is used to accomplish this 
task. 

Mass Balancing for Spinning Vehicle . When the cargo 
elevator is in motion, toe shitting of a large mass causes 
undesirable motions of the spinning vehicle. To balance 
this effect, an accelerometer detects the resulting 
accelerations, and a balance mass is moved to cor- 
rect both the center of mass shift and the cross pro- 
ducts of inertia of the vehicle. 

Miscellaneous Command Subroutines . The same sub- 
routine used to command the undeformed length of the 
Connecting Structure is also used to command the posi- 
tion (along each of the three axes) of every moving 
point mass except the balance mass (which is governed 
by toe balance mass control system). There are a 
maximum of seven such moving masses which may be 
used to simulate the cargo elevator, crew members, 
etc. Another subroutine is used to pre specify the 
fluid velocity in the piping system, i.e. , an idealized 
pump is assumed. Automatic shutoff occurs when a 
reservoir is either emptied or filled. This subroutine 
also generates the spin-rate command used as a 
reference input by the spin-rate control system. 




V, Miscellaneous Program Capabilties 

For reference, the program computes the position 
of the Space Station center of mass in inertial coor- 
dinates, the total angular momentum vector projected 
onto inertial coordinates, and the total system kinetic 
energy. The output includes ( 6 } and ( <t> * } , the linear and 
angular displacements, respectively, of the Connecting 
Structure attachment point on the Laboratory relative 
to the attachment point on the Counterweight, in addi- 
tion, the program can compute the internal resultant 
force and torque vectors on any surface that separates 
the structure into two free bodies. This is accomplished 
by the acceleration method, which appears to be more 
accurate than the stiffness-damping matrix method 
for the case where a truncated number of modes repre- 
sents the deformation of the vehicle. 

VI. Program Checkout 

Several teat problems were run to check the pro- 
gram. In the test configuration, the Laboratory was 
Idealized using eight structural masses and the Coun- 
terweight idealization had five structural masses. A 
telescopic boom was used as the Connecting Structure. 
First the entire vehicle was commanded to be rigid, 
and the time history was found to bn identical to results 
obtained using the well-known analytical solutions for 
a rotating rigid body. Next, a separate program was 
written by an independent programmer using the theo- 
retical expressions of Reference 1. However, this 
check program was much simpler than the main pro- 
gram since the expressions were applied only to the 
test configuration. Also, the time- and storage-saving 
manipulations used in the main program were not used 
In the test program; Several test runs were made, in- 
cluding vib lotion excited by: giving each variable and 
Its derivative a different initial condition while rotating; 
fluid motion; point masses in motion; and a retraction 
maneuver. Time histories of the two programs agreed 
in each case. 

The preprogram which synthesizes the modes of the 
Space Station was also tested by comparing results for 
several configurations with independently obtained 
results. The independent results were obtained by com- 
puting a mass and stiffness matrix for the entire struc- 
ture by a standard technique and applying an existing 
eigenvalue program. Agreement was achieved in each 
case. 

VII. Numerical Results Demonstrating 
the Computer Programs 

Configuration. The NASA Langley Research Center 
provided Grumman with the configuration shown in 
Figure 2 for the purpose of demonstrating the com- 
puter program. The Laboratory is composed of a 
central core module with two relatively rigid appen- 
dages and two solar panels which arc very flexible. 

The Connecting Structure is a telescopic beam which 
can be fully retracted. The Counterweight is com- 
posed of three relatively rigid modules. The mass 
of the entire vehicle is approximately 77,600 kg, and 
its deployed overall length is approximately 78.4 m. 

For additional descriptive details of this configuration 
see Reference 1, 

First, the modes of the core and appended modules 
were computed by Grumman, based upon stiffness 
and mass properties generated by North American 
Rockwell. The modes of the solar panels were 
computed by Fairchild Industries and Wolf Research. ' 
These modes were used as input data to the pre- 


program in order to obtain the free-free modes of the 
Laboratory and the C ounterweight. The Counterweight 
was relatively rigid. The lowest flexible Counterweight 
frequency was 6. 851 Ilz, whereas the sixth Laboratory 
frequency was .382 Hz. It was therefore decided to 
idealize the Counterweight as a rigid body in the time- 
history computer program. The Laboratory was 
idealized using 72 mass points. Six flexible laboratory 
modes were then used to reduce the number of 
numerical-integration coordinates. Because of the 
relatively high flexibility of the solar panels, most of 
the motion in these modes is solar- panel motion. 

A total of 18 coordinates are used in the time-history 
runs. These are the six laboratory modes, six 
rigid -body coordinates locating the mean axes for the 
Laboratory, and six rigid-body coordinates for the 
Counterw eight. 

The connecting structure was assumed to be an 
elastic tubular beam with uniform characteristics 
per unit length. To approximate telescopic charac- 
teristics of the beam during deployment and re- 
traction, new stiffness properties are computed at 
each time interval based on the beam length. 

All control system jets, the CMG, and the control- 
system sensor package were mounted on the Labora- 
tory core module. 

A detailed description of the configuration is 
presented in Section 6. 0 of Reference X. 

Selection of Runs. The runs which will be discussed 
are attitude control, deployment, spin-up, wobble 
control, elevator motion with mass balancing, and 
fluid being pumped between reservoirs. It is also 
possible to use the program to study many of these 
effects occurring simultaneously; however, this was 
not done in the present study. Since the runs 
performed were selected primarily to demonstrate 
the capability of the computer program, there was 
no attempt to optimize any of the operation param- 
eters such as control-system gains. Also, in the 
case of the deployment and spin-up maneuvers, the 
jet thrusts were increased to unrealistic values in 
order to complete the maneuvers within 45 sec, 
thereby saving computer time. 

Although the runs described herein arc primarily 
flexible-vehicle runs, in each case the rigid-vehicle 
run was also performed. Rigidization was accom- 
plished by using the program’s constraint option. 

Attitude Control. During the attitude-control 
maneuver, the system was not rotating and Uic 
Connecting Structure was fully retracted. Tho 
attitude control system developed for the rigid 
vehicle required no modification for use on the 
flexible-vehicle idealization. The three components 
of the Euler angles {y\ orienting the vehicle are 
shown in Figure 4. Curves for a rigid-body run 
were overlayed with the flexible body curves of 
Figure 4 and no difference could be discerned. 

The Space Staff on was initially tilted so that each 
Euler angle was .01745 rad (1.0 deg). The control 
system then reduces each angle to the commanded 
value of zero. Similar behavior occurs along each of 
the three axes. The jets first apply a torque to begin 
correcting the attitude angle. Then a torque is 
applied in the negative direction to slow down the 
Space Station’s angular rate. Linear deformation 
in the Connecting Structure is shown in Figure 5. 

Since the Connecting Structure is fully retracted, 
the illustrated deformation actually respesents the 
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Fig. 5 Cantilever Deflection | b 1 of Connecting Structure 
During Attitude Control Maneuver 


very small deflections of the relatively stiff Con- 
necting Structure docking hatch. This connection 
is only approximately represented as discussed in 
Section 4.4. 2.1 of Reference 1. 


Deployment with Spin-Rate Hold. Figure G illus- 
trates a deployment maneuver at .2 rpm. All motion 
in this run is in the spin plane. Figure GA shows the 
prespecified motion of the undeformed Connecting 
Structure. The Connecting Structure deploys from 
an undeformed length t of 0 to 42.822 m. If the 
spin control system were not operational, the spin 
rate would decrease to maintain constant momentum; 
however, the command to maintain a constant spin 
rate was given during this run. To accomplish the 
deployment within 40 sec, the jet thrust was increased 
from a nominal value of 222. 4 N (50 lb) to 4448. 2 N 
(1000 lb). Figure 6 B shews the axial deformation in 
the beam during deployment. Throughout the maneu- 
ver, the Blow spin rate causes a slight centrifugal 
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C CANTILEVER SENDING DISPLACEMENT OF CONNECTING 
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Fig. 6 Deployment Maneuver at .2 RPM 


force. Therefore, there is a small tendency for the 
end bodies to separate naturally; however, this would 
occur at a very slow rate. The deployment command 
is given at t ^ 1 see, and the Laboratory and the 
Counterweight are pushed apart (£<>3 >0 a 9 shown in 
Figure GA) causing compression in the beam. When 
the maximum velocity £ 03 = 1.270 m/sec (50 in. /sec) 
is reached, and deployment proceeds at constant 
velocity, the beam is expanded slightly by the centrif- 
ugal force. At approximately 35 sec, the deceleration 
begins ( £ q;j < 0 ) and the expansion in the beam is 
increased significantly. Deployment ends at about 
40 see. The final expansion is much larger than 
the initial compression, mainly because the Con- 
necting Structure Is more flexible when more of it 
is deployed. For the same reason, the transient 
vibration occurs at a lower frequency when deploy- 
ment is near completion. Figure 6 C illustrates 
the bending of the Connecting Structure during deploy- 
ment. This bending occurs partially aa a result of 
the Coriolis forces, but primarily it is due to the spin 
jets torquing the Laboratory to maintain a constant 
angular velocity. 


Spin-Up. Figures 7 and 8 illustrate a spin-up ' 
maneuver. The Space Station is initially rotating at 
.2 rpm and, at 5 sec, the command is given to in- 
crease the spin speed to 4 rpm. All motion in this 
run is in the spin plane. In order to accomplish the 
maneuver in 40 sec, the jets on the Laboratory were 
increased from 222.4 N (50 lb) to 66.7 kN (15,’000 lb). 
Figure 7 A illustrates the increase of the spin speed, 
and Figure 7B illustrates the corresponding increase 
in the axial extension of the Connecting Structure due 
primarily to the centrifugal force. Figure 7C shows 
the bending in the Connecting Structure during the 
spin-up maneuver. Unusually high bending deforma- 
tions occur as a result of the large torques on the 
Laboratory generated by the increased jet thrusts. 
Figure 8 shows the largest component of the deforma- 
tion at mass point 17, a very flexible point on a solar 
panel, and the largest component of the force exerted 
by the solar panel on the core module at the root of 
the panel (sec Figure 2). 

Quiescent State, When the Space Station is rotating 
in its nominal state of pure spin (i.e. , in a state of 
pure rotation about the Xj axiB with no vibration), 
constant elastic deformations occur due to the cen- 
trifugal force. This state Is known as the quiescent 
state. During the . runs which were made when the 





C. BENDING DISPLACEMENT OF CONNECTING STRUCTURE tN SPIN PLANF 


Fig* 7 Spin Speed and Cantilever Deflection of 

Connecting Structure During Spin-Up Maneuver 
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Fig. 8 Selected Deflection and Load Components 
During Spin-Up Maneuver 


Space Station was rotating at Us nominal spin speed 
of 4 rpm, the initial conditions are a variation from 
the quiescent state. Before making these runs, the 
quiescent state deformations were determined by 
setting all of the damping coefficients (for both die 
Laboratory and the Connecting Structure) to 80% of 
their critical values. A short run was made, and the 
deformations rapidly damped to their quiescent 
values. The quiescent deformations were highest 
at certain points on the solar panels. As an example, 
the deformation in the X 3 direction at mass point 17 
(qi? g) was approximately 540 mm (21.4 in.). 


Wobble Control , Figure 9 illustrates the perfor- 
mance of the wobble control system. ‘Initially, the 
deformations were set to their quiescent values, and 
the second component of (cd x > was given a wobble 
component of . 001 rad/sec. Up to approximately 
27 sec the curves are essentially identical to a run 
made for a rigid Space Station. This indicates the 
usefulness of the mean axes; one reason that they 
were used was that they move at the average motion 
of the deformed system. After 27 see, some small 
higher -frequency oscillations predominate due to 
elastic vibration. The only modification required to 



Fig. 9 Components of Angular Velocity |w X ) Illustrating 
Wobble Control System 


the CMG control system was to increase the amount 
of wobble at which the CMG stops trying to control 
wobble. Before this modification was made, the CMG 
sensor reacted to the residual vibration and the CMG 
continued to operate hi the wobble-damping mode 
throughout the entire run. 


Elevator Motion with Balance Mass Control . In 
the runs described in this paragraph, a 4,530 kg 
(10,000 lb) elevator and a 2,270 kg (5, 000 lb) balance 
mass are initially located on the X 3 axis, at X 3 -19.4 
and H.9 m, respectively. The elevator moves towards 
the balance mass at t=5 seconds. All rigid-body and 
flexible motion is in the spin plane. Initially, the 
Space Station is rotating in the quiescent state. The 
ffrst balance-mass control system which was devel- 
oped operated properly on a rigid-vehicle idealiza- 
tion; however, using the current computer program 
it was found that vehicle flexibility coupled undesir- 
ably with the control system. During this run the 
vehicle vibrated at a high-frequency vibration 
associated with the Connecting Structure axial mode, 
and the overall response of the balance mass waB 
sluggish and unstable. The satisfactory control, 
illustrated in Figure 10 was achieved after 
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Fig. 10 Performance of Balance Mass Control System 



considerable modification of the control law. The 
control system Is required to balance an elevator 
motion of 7.00 m. The curve showing the error fn the 
position of the Space Station center of mass indicates 
that there is a lag in the response of the balance mass; 
however, by the end of the run the Space Station is 
balanced. The deformations in the Connecting 
Structure during this run are shown in Figure 11. 

Note that there is an initial vibration although the 
elevator does not begin its motion until fi sec. This 
initial vibration results because the initial quiescent 
deformations were obtained from a run where the 
elevator and balance masses were not present. 

The main Connecting-Structure bending effects are 
caused by the Coriolis forces exerted by the moving 
masses on the Laboratory and by the spin jets 
which torque the Laboratory to maintain the com- 
manded constant spin speed. The deformation at 
mass point 17, which is on a solar panel, and the 
internal force exerted by the solar panel at its root 
are shown in Figures 12A and B . 


ill Figure 13. Pumping proceeds until this reservoir is 
empty at t = 35.7 sec, when the pump suddenly shuts 
down. After shutdown, fluid remains in the pipeline. 

Alt motion during this run occurs in the spin plape. 
Figure 14 shows the deformations of the Connecting 
Structure. The bending which occurs during pumping 
is illustrated in Figure 15. The primary reason for 
this bending is that the resultant of the Coriolis forces’ 
exerted by the fluid acts to the left of the cm of the 
Laboratory as illustrated in the figure. This force also 
tends to slow the spin speed of the Laboratory slightly 
(from .41H9 to .4149 rad/sec.). The control system 
is not present during this run; therefore, the spin 
jets do not turn on to correct the spin speed. 
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Fig. 11 Cantilever Deflection of Connecting Structure 
During Mass Balancing 
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Fig. 14 Cantilever Deflection of Connecting Structure 
During Fluid Motion 
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Fig. 12 Selected Deflection and Load Components 
During Mass Balancing 

Fluid Pumped Through the Laboratory. In this 
run only, the Laboratory contained a fluid system con- 
sisting of a pipeline connecting two reservoirs. Unlike 
the system illustrated in Figure 2, the centerlines of 
both reservoirs and the pipe were located on the 
axis. Pumping begins at 2 Bee. The fluid velocity and 
height of the fluid In the .emptying reservoir are shown 
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Fig. 15 Sketch Illustrating Bending (Exaggerated) in 
Connecting Structure During Fluid Motion 
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